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Abstract 

We present a new measurement of the inclusive and differential production cross sections of J/V' 
mesons and 6-liadrons in proton-antiproton collisions at i/s = 1960 GeV. The data correspond to 
an integrated luminosity of 39.7 pb^^ collected by the CDF Run II detector. We find the integrated 
cross section for inclusive J/V' production for all transverse momenta from to 20GeV/c in the 
rapidity range \y\ < 0.6 to be 4.08 ± 0.02(stat)]to 33(^2/5*) yub. We separate the fraction of J/tp 
events from the decay of the long-lived 6-hadrons using the lifetime distribution in all events with 
ppiJ/ip) > 1.25GeV/c. We find the total cross section for 6-hadrons, including both hadrons and 
anti-hadrons, decaying to J/ip with transverse momenta greater than 1.25GeV/c in the rapidity 
range \y{J/ip)\ < 0.6, is 0.330 ± Om5{stat)to^ll{syst) jih. Using a Monte Carlo simulation of the 
decay kinematics of 6-hadrons to all final states containing a J/i/j, we extract the first measurement 
of the total single 6-hadron cross section down to zero transverse momentum at s/s = 1960 GeV. 
We find the total single 6-hadron cross section integrated over all transverse momenta for 6-hadrons 
in the rapidity range \y\ < 0.6 to be 17.6 ± 0.4:(stat)~^2 3{^y^i) M^- 



PACS numbers: 13.85.Qk, 13.20.Gd, 14.40.Gx, 12.38.Qk 



I. INTRODUCTION 



The production of both charmonium mesons and bottom-flavored hadrons (referred to as 
6-hadrons or Hh in this paper) in proton- antiproton coUiders has sustained continued interest 
over the last several years. There are three major sources of the J/t/j mesons: directly 
produced J/ip, prompt decays of heavier charmonium states such as ^Pi state Xci and 
state Xc2, and decays of 6-hadrons. Early hadroproduction models of quarkonium states 
could not describe the cross section of directly produced J/ -0 mesons. These models under- 
predicted the measurements by a factor of approximately 50, and did not adequately describe 
the cross-section shape [1] . With the advent of the effective field theory, nonrelativistic QCD 
(NRQCD) [2], better theoretical descriptions of quarkonium production became possible. 
Within the NRQCD factorization formalism, the color-octet model provides a means to 
bring theory into better agreement with data [3, 4]. The fundamental idea of this model 
is that while a (cc) meson has to be in a color-singlet state, the initially produced quark- 
antiquark pair does not. One can produce, for example, a (cc) pair in a color-octet state 
which can then produce a color singlet ^Si J/ip meson by single-gluon emission. This is 
done at the cost of adding a small number of parameters to the theory that currently must 
be determined by experiment. While the color-octet model can accommodate a large cross 
section, strictly speaking it does not predict it. There are other deficiencies of the NRQCD 
formalism; for example, NRQCD expects that the spin alignment to be predominantly in 
the transverse state for the prompt J/ ip mesons with large transverse momenta, a prediction 
that is not in agreement with subsequent measurement [5]. 

Previous prompt, direct, and inclusive J/ip cross-section measurements [1] from CDF 
required a minimum transverse momentum of 5GeV/c on the J/ip although greater than 
90% of the cross section has been expected to lie below this point. In this paper we present 
the first measurement of the inclusive central J/ ip cross section over a much larger range of 
transverse momenta from zero to 20GeV/c. The J/ifj mesons are reconstructed from the 
decay channel J/ip ^ A*^/^~- The measurement was made possible by improving the CDF 



di-muon trigger capability to be sensitive to J/V' with zero transverse momenta. 

A significant fraction of J/%1) mesons produced at the Tevatron come from the decays of 
6-hadrons. In this experiment, we use the large sample of if^ — > J/%l)X events to measure 
the inclusive 6-hadron cross section. The previous Tevatron measurements [6-12] of the 
6-hadron cross section in proton- antiproton collisions at = 1800 GeV were substantially 
larger (by a factor of two to three) than that predicted by next-to-leading order (NLO) 
QCD calculations [13-15]. This was particularly puzzling since the UAl measurements at 
^/s — 630 GeV [16] did not show such a marked departure from the NLO QCD calculations. 
Several theoretical explanations were suggested: higher-order corrections are large, intrinsic 
kr effects are large [17], extreme values of the renormalization scales are needed, or new 
methods of resummation and fragmentation are required [18-20]. Theories of new and 
exotic sources of 6-hadrons have also been proposed [21]. Since the earher Tevatron results 
covered only 10-13% of the inclusive pr spectrum, it was not evident whether the excess was 
due to an overall increase in the 6-hadron production rate or a shift in the spectrum toward 
higher pt- 

An inclusive measurement of 6-hadron production over all transverse momenta can help 
resolve this problem. Bottom hadrons have long lifetimes, on the order of picoseconds [22] , 
which correspond to flight distances of several hundred microns at CDF. We use the measured 
distance between the J /ip decay point and the beamline to separate prompt production of 
charmonium from 6-hadron decays. The single 6-hadron cross section is then extracted 
from the measurement of the cross section of J/tl) mesons from long-lived 6-hadrons where 
single differentiates the cross section from the 6-hadron cross section referring to h and h 
hadrons which is a factor of two bigger. In this paper, we present the first measurement of 
the inclusive single 6-hadron cross section at ^/s — 1960 GeV measured over all transverse 
momenta in the rapidity range \y\ < 0.6. 



II. DESCRIPTION OF THE EXPERIMENT 



A. The Tevatron 

The Fermilab Tevatron is a 1 km radius superconducting synchrotron. Thirty-six bunches 
of 980 GeV protons and antiprotons counter-circulate in a single ring and collide at two in- 
teraction points (where the CDF and DO detectors are located) every 396 ns. The transverse 
profile of the interaction region can be approximately described by a circular Gaussian dis- 
tribution with a typical RMS width of 30 /im. The longitudinal profile is also approximately 
Gaussian with a typical RMS of 30 cm. 

For the data used in this analysis, instantaneous luminosities were in the range 0.5 to 
2.0 X 10^^ cm~^s~^. At these luminosities, typically there was only a single collision in a 
triggered event. 

B. The CDF Detector 

In the CDF detector [23, 24], a silicon vertex detector (SVX II) [25], located immediately 
outside the beam pipe, provides precise three-dimensional track reconstruction and is used 
to identify displaced vertices associated with b and c hadron decays. The momentum of 
charged particles is measured precisely in the central outer tracker (COT) [26], a multi- 
wire drift chamber that sits inside a 1.4 T superconducting solenoidal magnet. Outside the 
COT are electromagnetic and hadronic calorimeters arranged in a projcctive-tower geometry, 
covering the pseudo-rapidity region \r)\ < 3.5. Drift chambers and scintillator counters in the 
region jTyj < 1.5 provide muon identification outside the calorimeters. In the CDF coordinate 
system, 9 and are the polar and azimuthal angles, respectively, defined with respect to the 
proton beam direction, z. The pseudorapidity r] is defined as — lntan(6'/2). The transverse 
momentum of a particle is pr — psin(^). 

The portion of the sihcon detector systems used in this analysis is the SVX II detector. 
The SVX II consists of double-sided micro-strip sensors arranged in five concentric cylindrical 



shells with radii between 2.5 and 10.6 cm. The detector is divided into 3 contiguous five-layer 
sections along the beam direction for a total z coverage of 90 cm. Each barrel is divided into 
twelve azimuthal wedges of 30° each. Each of the five layers in a wedge is further divided 
into two electrically independent modules called ladders. There are a total of 360 ladders in 
the SVX II detector. The fraction of functioning ladders was increasing from 78% to 94% 
during the period between February 2002 and July 2002 in which the data used in this paper 
were taken while the SVX detector was being commissioned. 

The COT is the main tracking chamber in CDF. It is a cylindrical drift chamber seg- 
mented into eight concentric superlayers filled with a mixture of 50% Argon and 50% Ethane. 
The active volume covers \z\ < 155cm and 40 to 140cm in radius. Each supcrlayer is sec- 
tioned in 4> into separate cells. A cell is defined as one sense plane with two adjacent 
grounded field sheets. The sense plane is composed of 40 iim gold-plated tungsten wires, 
twelve of which are sense wires. In the middle of the sense planes, a mechanical spacer made 
of polyester /fiber glass rod is epoxied to each wire to limit the stepping of wires out of the 
plane due to electrostatic forces. The main body of the field sheets is 10 iim gold-coated 
mylar. The field sheets approximate true grounded wire planes much better than the arrays 
of wires which have often been used in wire chambers including the predecessor to the COT. 
Use of the field sheet also results in a smaller amount of material within the tracking volume, 
and allows the COT to operate at a much higher drift field than is possible with an array of 
wires. The eight superlayers of the COT alternate between stereo and axial, beginning with 
superlayer 1, which is a stereo layer. In an axial layer, the wires and field sheets are parallel 
to the z axis, and thus provide only r-cf) information. In stereo layers, the wires and field 
sheets are arranged with a stereo angle of ±2° and provide z information in addition to r-0. 

The CDF central muon detector (CMU) [27] is located around the outside of the central 
hadron calorimeter at a radius of 347 cm from the beam axis. The calorimeter is formed 
from 48 wedges, 24 on the east (positive z), and 24 on the west (negative z), each wedge 
covering 15° in 0. The calorimeter thickness is about 5.5 interaction lengths for hadron 
attenuation. The muon drift cells with seven wires parallel to the beamline are 226 cm long 



and cover 12.6° in (j). There is a 2.4° gap between drift cell arrays, giving a (j) coverage 
of 84%. The pseudorapidity coverage relative to the center of the beam-beam interaction 
volume is 0.03 < \rj\ < 0.63. Each wedge is further segmented azimuthally into three 4.2° 
modules. Each module consists of four layers of four rectangular drift cells. The sense wires 
in alternating layers are offset by 2 mm for ambiguity resolution. The smallest unit in the 
CMU, called a stack, covers about 1.2° and includes four drift cells, one from each layer. 
Adjacent pairs of stacks are combined together to form a two-stack unit called a tower. A 
track segment detected in these chambers is called a CMU stub. 

A second set of muon drift chambers is located behind an additional 60 cm of steel (3.3 
interaction lengths) . The chambers are 640 cm long and are arranged axially to form a box 
around the central detector. This system is called the CMP, and muons which register a 
stub in both the CMU and the CMP are called CMUP muons. 

Luminosity is measured using low-mass gaseous Cherenkov luminosity counters 
(CLC) [28, 29]. There are two CLC modules in the CDF detector installed at small an- 
gles in the proton and antiproton directions. Each module consists of 48 long, thin conical 
counters filled with isobutane gas and arranged in three concentric layers around the beam 
pipe. 

C. Muon Reconstruction 

The starting point for the selection of J/t/j — > candidates is the reconstruction 

of two oppositely charged muons. Muons are reconstructed from tracks measured in the 
tracking chambers matched to the stub positions in the muon detectors. 

1. Charged Particle Tracking 

Track reconstruction begins in the COT. The first step in the pattern recognition is the 
formation of line segments from hits in each superlayer. Line segments from the axial layers 
that are tangent to a common circle are linked together to form a track candidate and the 



hit positions are fit to a circle. Line segments in stereo layers are then linked to the 2- 
dimensional track and a helix fit is performed. The transverse momentum resolution of the 
COT is measured using cosmic ray events to be 

= 0.0017 [GeY/c]-\ (1) 

Pt 

The next step is to extrapolate each COT track into the SVX II and add hits that are 
consistent with the track. A window around the track is established based on the errors on 
the COT track parameters. If a hit in the outer SVX II layer lies within the window, it 
is added to the track. A new track fit is then performed, resulting in a new error matrix 
and a new window. This window is then used to add hits from the next SVX II layer, and 
the procedure is repeated over all layers. If no hit is found within the search window, the 
algorithm proceeds to the next layer. There may be multiple track candidates with different 
combinations of SVX II hits associated with one COT track. In this case, the track with 
the largest number of SVX II layers with hits is chosen. A COT-SVX II track is formed 
only if at least three r-0 hits in the SVX II are associated with the original COT track. An 
averaged impact parameter resolution of 34 /xm is achieved using hit information measured 
in SVX II for muon tracks with pt around 1.5 GeV/c. 

2. Muon Identification 

In the first stage of muon identification, the measured drift times for hits in the muon 
chamber drift cells are converted to drift distances. Hits in alternate layers that are within 
7.5 cm of each other are used to form linear track segments. This distance corresponds to a 
maximum angle relative to the radial direction in the chamber of 65°. The remaining pair of 
layers is then searched for hits within 0.5 cm of the line segment. The procedure is iterated 
and the optimal set of hits is found. The segment resulting from a least-square fit to these 
hits is called a "stub" . Hits are required in at least 3 of the 4 layers to form a stub. 

Stubs reconstructed in the CMU are matched to tracks with a minimal pt of 1.3 GeV/c. 
The tracks are extrapolated to the CMU after using a simplified geometry model to track 



the muon candidate's motion in the non-uniform magnetic field of the calorimeter. The 
distance, Ar0, in the r-0 plane between the track projected to the muon chambers and the 
muon stub is required to be less than 30 cm. The track is required to point to the same side 
(east or west) of the detector that the stub is in unless the muon candidate track is within 
20cm of the center of the detector. 

D. Triggers 

CDF uses a three-level trigger system [24]. At Level 1, data from every beam crossing is 
stored in a pipeline capable of buffering data from 42 beam-crossings. The Level 1 trigger 
either rejects the event or copies the data into one of four Level 2 buffers. During the data- 
taking period for this analysis, the global Level 1 accept rate was approximately 10 kHz 
corresponding to a rate reduction factor of approximately 170. 

At Level 2, a substantial fraction of the event data is available for analysis by the trigger 
processors which require approximately 25 /xs per event. During the period the data for this 
analysis were taken, the L2 accept rate was approximately 200 Hz, for a rejection factor of 
approximately 50. 

Events that pass the Level 1 and Level 2 selection criteria are then sent to the Level 3 
trigger [30], a cluster of computers running a speed-optimized reconstruction code. Events 
selected by Level 3 are written to permanent mass storage. During the period the data 
for this analysis were taken, the global Level 3 accept rate was approximately 40 Hz, for a 
rejection factor of approximately 5. 

For the cross-section measurement, we require events with two muon candidates identified 
by the Level 1 trigger. In Level 1, track reconstruction is done by the eXtremely Fast Tracker 
(XFT) [31]. The XFT examines COT hits from the four axial superlayers and provides r-0 
tracking information. The line segments are identified in each superlayer and linked using 
predetermined patterns. The XFT requires that each fine segment contains hits found on at 
least ten of a possible twelve anode wires in each axial superlayer. The XFT finds tracks with 



Pt > 1.5GeV/c. It subdivides the COT into azimuthal sections of 1.25° each and places a 
track into a given section based on its 4> position at superlayer 6 (r = 105.575 cm). If more 
than one track candidate is found within a given section, the XFT return the track with the 
highest Pt. The XFT passes the tracks it finds to the eXTRaPolation unit (XTRP). The 
XTRP extrapolates an XFT track's trajectory to the CMU where a stub should be found if 
it is a muon, taking into account the path of the track in the magnetic field and the multiple 
scattering of muon in the calorimeter. The XTRP then passes the search window to the 
muon trigger crate, which looks for CMU stubs within the search window. A Level 1 CMU 
stub requires that there be hits on both even layers or both odd layers of one 1.05° stack 
of the CMU with a drift time difference At less than 396 ns. The twelve stacks in each 15° 
wedge of the CMU are mapped in pairs to six trigger towers to match the granularity of the 
XTRP extrapolation. If a muon stub is found within the search window, it is considered a 
Level 1 muon. In order to fire the di-muon trigger, two muon candidates must be found, 
separated by at least two CMU trigger towers. There is no requirement that the muons have 
opposite charge at Level 1. During the data taking period in which the di-muon sample used 
for this analysis was obtained, there is no additional selection imposed on muons at Level 2 
and event is passed to Level 3 directly from Level 1. 

At Level 3, the muons are required to have opposite charge, and to have an invariant 
mass between 2.7 and 4.0GeV/c^. In addition, both muon tracks are required to be within 
5 cm in zq, where zq is the z coordinate of the muon track at its distance of the closest 
approach in the r-0 plane to the beam axis. For a portion of the data sample considered in 
this analysis, there is a requirement that the opening angle in r-0 between the di-muons be 
less than 130°. 

E. Luminosity 

The CLC counters monitor the average number of inelastic pp interactions in each bunch 
crossing. The inelastic pp cross section has been measured to be ~ 60 mb by several 



experiments at ^/s = 1800 GeV [32-34]. The inelastic pp cross section at ^/s = 1960 GeV is 
scaled from previous measurements using the calculations in [35] . The rate of inelastic pp 
interactions is given by 

• /bC = CTin ■ (2) 

where C is the instantaneous luminosity, /i is the average number of inelastic pp interactions 
per bunch crossing, and /bc is the rate of bunch crossings. In this paper, we use data from 
the beginning of the CDF Run II operation where the average instantaneous luminosities 
were relatively low. 

The number of pp interactions in a bunch crossing follows Poisson statistics where the 
probability of empty crossings is given by Vq{^jl) — e~'^. An empty crossing is observed 
when there are fewer than two counters with signals above threshold in either module of the 
CLC. The measured fraction of empty bunch crossings is corrected for the CLC acceptance 
and the value of /i is calculated. The measured value of /i is combined with the inelastic 
pp cross section to determine the instantaneous luminosity using Equation 2. Because this 
method depends only weakly on the CLC thresholds, it functions particularly well at low 
luminosities where the probability of empty bunch crossings is large. The systematic error 
on the luminosity measurement is estimated to be 6%. 

In CDF Run II, only runs with greater than 10 nb~^ integrated luminosity are considered 
for analysis. Runs with good operating conditions in the detector are tagged by the online 
shift crews. Data from those runs are examined to exclude ones with COT, muon or trigger 
hardware problems. For the measurement presented in this paper, the data collected from 
February to July 2002 was used. This sample corresponds to a total integrated luminosity 
of 39.7 ±2.3 pb-^ 

For J/i/j candidates with transverse momenta in the range to 2 GeV/c, we use 14. 8 ±0.9 
pb~^ of our data sample, which corresponds to that fraction of the data collected when no 
cut on the di-muon opening angle in the Level 3 trigger was used. 



III. DATA SELECTION AND RECONSTRUCTION 



A. Data selection 

The events selected by the J/tp trigger are reconstructed offline, taking advantage of the 
most refined constants and algorithms. We reconstruct J/%1) ^ A*^A*~ decays by selecting 
events with two oppositely charged muon candidates reconstructed in the COT and CMU 
detectors. The J/i]j — > sample used for this analysis was collected using the CMU 

di-muon triggers. Events are required to have satisfied the Level 1 and Level 3 di-muon 
trigger criteria. 

In addition to the default muon selection criteria outlined earlier, we require a px inde- 
pendent track-stub matching criterion x^(Ar0) < 9. A track-stub matching quality criterion 
X^(Ar0) with a one degree of freedom is calculated from Arcf) and the expected multiple 
scattering for a track of given Pt obtained from a GEANT simulation [37] of the CDF Run 
II detector material. We require both muons to have transverse momenta px > 1.5 GeV/c 
as measured offline. The trigger requirements are verified for the offline-reconstructed can- 
didates. In addition, each CMU stub matched to a triggered stub must lie within the XTRP 
search window set by the Level 1 triggered track. Furthermore, track momentum is cor- 
rected for energy loss due to specific ionization and multiple scattering according to our 
accounting of the detector materials. We calculate the J/ ip candidate invariant mass from 
the four-momenta of the two muons. For a portion of the data sample under consideration, 
a temporary hardware problem with the di-muon logic caused the trigger to exclude J/ip 
events where both muon stubs fell in the range of 240 — 270°. Therefore, we exclude 
J/ip events where both muons fall in that (p region and account for this in the detector 
acceptance. We also reject J ftp candidates if one of the tracks passes within 1.5 cm of the 
center of any COT wire planes, where the trigger efficiency is difficult to model because of 
the distortion of the electric field due to the mechanical spacers. This exclusion is accounted 
for in the acceptance calculation. The muon reconstruction efficiency is measured in each of 
the 48 CMU detector wedges. We find that the hit efficiency in the CMU wedge on the west 



side of the detector covering the region 240° < < 255° is lower due to a known hardware 
problem and exclude J/ ip events where either muon stub is reconstructed in this wedge. As 
shown in Fig. 1, there are 299800 ± 800 J/tp events that passed these selection conditions. 
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FIG. 1: Mass distribution of reconstructed di-muon J/ip candidates. The points are data. The solid 
line is the fit to the signal approximated as a double Gaussian and a linear fit for the background. 
The hatched region is the fitted background. The fit gives a signal of 299800 it 800 J/tp events 
with an averaged mass of 3.09391 it 0.00008 GeV/c^ obtained and an average width of 0.020 it 
0.001 GeV/c^ mainly due to detector resolution. The uncertainties here are statistical only. 

To determine the yield in each J/ip px bin, the di-muon invariant mass distribu- 
tions are fitted using invariant mass line shapes including the radiative tail from internal 
bremsstrahlung obtained from a tuned hit-level COT simulation. The simulated J/ ip events 
are decayed using the J/ip radiative decay model in the QQ decay package [38]. The COT 
hit multiplicity per track is tuned to match the data as closely as possible. The COT hit 
resolution is then tuned to find the best ^ binned fit to the data using the Monte 

Carlo invariant mass line shape for the signal and a polynomial shape for the background. 
Finally, energy loss and multiple scattering in material encountered before the COT are 
modeled. The energy loss in the silicon material is scaled until the peaks of the di-muon 



invariant mass distribution in different pt ranges in data and from tlie simulation match. 
Tlie order of the background polynomial used varies with the background shape in each 
J/%1) Pt range. A third-order polynomial is used for the momentum range 0-0.25 GeV/c, a 
second-order polynomial is used for the range 0.25-2.25 GeV/c, and a first-order polynomial 
(linear background) for transverse momenta greater than 2.25GeV/c. The fits to the in- 
variant mass distributions in four J/%1) px ranges are shown in Figs. 2-5. The J/%1) yields 
and the statistical uncertainties obtained from the fits in each pj- range are listed in the first 
column of Table I. The mass fitting qualities over the whole px bins are good as indicated 
from the fit probability shown in these Figures. We also examined the differences between 
counting the event numbers in the J/ iIj signal region ( 3.02 3.15 GeV/c^) to that predicted 
from the fitting functions of signal and background. The differences ranging from -|-9% in 
the lowest px bin to -1.3% in the high px bin are used very conservatively as the systematic 
uncertainties from the mass fitting. 

IV. ACCEPTANCE AND EFFICIENCY 
A. Monte Carlo Description 

We use the GEANT [37] Monte Carlo simulation software to estimate the geometric 
and kinematic acceptances. The variation of detector conditions in the simulation is set to 
match the data. J/%Ij events are generated starting with a kinematic distribution that is flat 
in rapidity and with a px distribution selected to best match the reconstructed data. The 
events are fully simulated. After the differential cross section is measured, we iterate and 
recalculate the acceptance and the central value of the cross section using the measured px 
distribution. The GEANT simulation is vahdated by comparing the resulting distributions 
of various kinematic quantities such as r^, px, the track-stub matching distance, and the z 
vertex distribution in reconstructed data and reconstructed Monte Carlo events. Differences 
in the data and Monte Carlo distributions are used to estimate the systematic uncertainties 
on the modeling of the CDF detector geometry in the simulation. 



TABLE I: Summary of the inclusive J/ij; cross-section analysis components. The values of the 
yield and statistical uncertainty from the fits are listed in the 2nd column. The acceptance values 
and the combined systematic and statistical uncertainties on the acceptance are listed in the 3rd 
column. In the 4th and 5th columns the trigger and track-stub matching efficiencies obtained from 
the mean of the distribution in each bin and the corresponding systematic uncertainties are listed. 
The sixth column lists the integrated luminosity used for each measurement. 
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FIG. 2: Invariant mass distribution of reconstructed J/ij: ^ ^u// events in the range pTip-fJ^) < 
0.25 GeV/c. The points with error bars are data. The sohd Une is the fit to the signal shape from 
the simulation and a third order polynomial for the background. The shaded histogram is the 
fitted background shape. The number of signal events and the fit probability of the binned 
fitting are also provided. 

B. Acceptance 

We correct the observed number of J/ip events for the detector acceptance and efficiency. 
The CMU muon detector covers the pseudo-rapidity range of | 77 |< 0.6. In this region the 
coverage of the COT is complete and the CDF detector acceptance is driven by the muon 
detector geometry and kinematic reach. The calorimeter acts as an absorber for the CMU 
detector which is therefore sensitive only to muons withpj- > 1.35 GeV/c. The arrangement 
of the four sense wires within the CMU chambers allows a lower bound on the transverse 
momentum of the muon to be calculated from the difference in drift times in sense wires on 
alternating layers. The At < 396 ns timing window is selected to be fully efficient for muons 
with pt > 1.5GeV/c. 

The acceptance is modeled as a function of both the reconstructed Pt{J/'4^) and rapid- 
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FIG. 3: Invariant mass distribution of reconstructed J/V' fJ-^i events in the range 1.25 < pxip-fj) < 
1.5GeV/c. The points with error bars are data. The sohd hne is the fit to the signal shape from 
the simulation and a second order polynomial for the background. The shaded histogram is the 
fitted background shape. 

ity y{J/ilj) and is defined as the ratio between the number of generated events A^^e" ^nd 
recontructed events N^'^'^, 



A{pt, y) 



iV--(pT(W),|l/(J/V^)|<0.6) 



(3) 



iV9-(p^(j/^),|^'(j/^)| <0.6)' 

where PTiJ/i^) and y' are the generated true values of the J/ip momentum and rapidity. 
The acceptance as a function of pt and y is shown in Fig. 6. 

The acceptance increases rapidly from 0.7% at px = 0.25GeV/c to 10% at 5GeV/c 
and 25% at 20GeV/c. The acceptance in the range 0.0-0.25 GeV/c is rapidly varying as 
a function of Pt{J/'4') aiid increases with decreasing momenta from 0.7% at Pt{J/4') = 
0.25GeV/c to 4% for J/tp mesons almost at rest {px < 50 MeV/c). The muon transverse 
momentum is required to be greater than or equal to 1.5GeV/c, which is close to one-half 
of the J/ ip mass, therefore when the J/ ip is at rest both muons are likely to be above the pt 
threshold. As soon as the J/ip receives a small boost, the probability is greater that at least 
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FIG. 4: Invariant mass distribution of reconstructed J/'0 — > /^A* events in the range 5.0 < pxip-fj) < 
5.5GeV/c. The points with error bars are data. The sohd line is the fit to the signal shape from 
the simulation and a linear background. The shaded histogram is the fitted background shape. 

one muon will be below the px acceptance threshold and the acceptance starts to decrease 
until the J/^p transverse momentum exceeds 0.25GeV/c. 

There is a small but non-zero acceptance at \y\ = 0.6 due to detector resolution and 
the size of the interaction region. J/ip Monte Carlo events generated with a fiat rapidity 
distribution in the range \y\gen < 1-0 and a pt distribution as described in Section IVA are 
simulated. The relative acceptance of events generated with \y\gen > 0.6 and reconstructed 
with \y\reco < 0.6, A', is calculated thus: 



N'^^{\y 


gon > 0.6, 


y 


roc < 0.6) 


N9^''{\y 


gen < 0.6) 



(4) 



where N^'^^{\y\^e^ > 0.6, \y\rcc < 0.6) is the number of J /ip events in the Monte Carlo sample 
with reconstructed |?/|reco < 0.6 and generated \y\gen > 0.6 and N^'^^{\y\gcn < 0.6) is the 
total number of events generated with \y\ < 0.6. The value of A' is found to be very small: 
A' = 0.00071 ± 0.00006(stat). A correction factor of (1 - A') = 99.93% is applied to the 
J/ip yield calculated in each Pt{J/4') bin. 
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FIG. 5: Invariant mass distribution of reconstructed J/V' — > fJ-^i events in the range 12.0 < pxip-fj) < 
14.0 GeV/c. The points with error bars are data. The sohd line is the fit to the signal shape from 
the simulation and a linear background. The shaded histogram is the fitted background shape. 

A 2-dimensional acceptance function was used for an event-by-event correction during 
the cross section calculation process. In Table I, the averaged acceptance values and the 
combined statistical and systematic uncertainties for each pt bin are given. Sources of 
systematic uncertainties studied are J/ip spin alignment, pt spectrum, CMU simulation and 
detector material description in GEANT simulation. 

Kinematic acceptance as a function of pt depends on the J/ip spin alignment. The 
normalized alignment distribution is given by 

3 



m 



1 + a cos 



(5) 



2(a + 3) 

where 6 is the angle between the muon in the J/ip rest frame and the direction of the 
J/ip in the lab frame [5] and a quantifies the spin alignment. The parameter a must lie 
in the range -1 to 1 and a = indicates no preferred spin alignment. The previous CDF 
measurements of the J ftp spin alignment parameter [5] are consistent with zero but could 
also be as large as 50% in some pt regions. The weighted mean of a measured in different 
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FIG. 6: Acceptance of J/^ /x/x events determined from a GEANT simulation of the CDF 
detector. The acceptance is shown as a function oipT{J/ip) and y{J/ijj). The acceptance as a 
function of priJ/i^) is measured integrated over |y| < 0.6 and the acceptance as a function of y is 
shown integrated over all pT- 



Pt ranges in [5] is used to determine the central value of the parameter a to be used for the 
acceptance. The value of o; = 0.13 ± 0.15 is used for the final acceptance values where the 



uncertainty is chosen to accommodate the variation in the previous CDF measurements and 
the extrapolation to pt = where a is expected to be zero. The uncertainty on acceptance 
due to spin ahgnment is largest in the lower momentum bins and decreases with increasing 
transverse momentum. We find the uncertainty is ~ 5% near pt — ^ and 2% in the region 
n <Pt< 20 GeV/c. 

To estimate the uncertainty from variations of the input transverse momentum spectrum, 
the acceptance is recalculated using a Monte Carlo sample generated using a flat distribution. 
The flat distribution is an extreme alternative from the nominal spectrum which is a fast 
falling function of pt- The fractional change in acceptance is taken as the uncertainty on 
the input transverse momentum distribution. The uncertainty is about 3% in the lowest 
momentum bin, less than 1% in the 0.25 to 3GeV/c bins, 1-2% in the 3 to 4GeV/c bins, 
and 2-4% in the 4 to 20 GeV/c bins. 

A systematic error of 1.0% from uncertainties related to the CMU chamber simulation is 
estimated by comparing event distributions in data and in Monte Carlo. The modeling of 
the CMU coverage in r-z plane, the wire efficiency differences between wedges in east and 
west and in different (f) sections, and beam position in z are found to be the major sources 
of the simulation uncertainties. 

There is a gap in CMU coverage in the central region of the detector in the r-z plane. 
The gap in coverage is approximately ±11 cm, measured at a radius of 347 cm. The fraction 
of muons falling in the gap region but still accepted by the CMU due to multiple scattering 
is compared between data and Monte Carlo. The deviation between the ratios in data and 
Monte Carlo is taken as the uncertainty in the modeling of the CMU fiducial volume in the 
center of the CDF detector. The uncertainty is found to be 0.20%. 

Several factors contribute to the difference in the numbers of J/%1) mesons with decay 
vertex in the opposite halves of the detector along z. These include the shift in the average 
primary vertex location towards positive z (east), the exclusion of the low efficiency wedge 
on the west side of the detector, and the uncertainty in the modeling of the z extent of 
the CMU detector, as well as the differences in the east and west chambers. We found 



a difference of 0.80% between data and Monte Carlo on tfie east-west asymmetry in tfie 
number of reconstructed events. 

The acceptance of the CMU detector obtained from the GEANT simulation does not 
include the differences in gain and efficiencies between wedges. The number of events recon- 
structed in each wedge in data and Monte Carlo is examined and the total number of events 
in Monte Carlo is normahzed to match the data. The standard deviation of the difference 
between the number of events reconstructed in each wedge between data and Monte Carlo 
is taken as the uncertainty on the CMU (j) acceptance. We find an uncertainty of 0.55% due 
to this source. 

Muons from J/%1) are required to have the Zq position to be within 90 cm of the center 
of the detector, |^o| < 90 cm. There is a small disagreement between the data and the 
Monte Carlo in the Zq{^) distributions due to inadequate modeling of the interaction region. 
This contribution to the systematic error is estimated from the difference between the ratios 
of data and Monte Carlo tracks with |^o| < 90 cm compared to all muons. We find an 
uncertainty of 0.28%. 

The material description of the CDF detector in GEANT determines the amount of 
energy loss from a muon track when it travels through the detector. Inside the tracking 
volume, the material description of the new silicon detector has the biggest impact on muon 
tracks in the low momentum range which is of special interest to this analysis. To estimate 
the systematic error on the acceptance from uncertainty of the detector material description, 
the SVX II material used in the simulation was varied by 10 % to 20%. The systematic 
uncertainty is taken as the difference between the acceptance values measured with different 
material scale factors and the nominal. The uncertainty is largest in the low momentum 
bins where it is around 5%. 

The systematic uncertainties on acceptance are summarized in Table IV. The size of 
the uncertainties from J/ if) spin alignment, J/ if) pr spectrum and detector material descrip- 
tion depends on the muon px range as expected while the uncertainty from muon detector 
simulation is same for all pr ranges of interests in the analysis. 



C. Data Quality 



The yield, mean, and resolution of the J /ip invariant mass peak were monitored over the 
period of the data taking to evaluate the detector performance. The number of J/ if) mesons 
reconstructed is normalized by the integrated luminosity of each run. We identify outlying 
runs which may have additional hardware or trigger problems that have been undetected by 
the standard offline validation procedures. Runs with J/ if) yields different by 4a" from the 
average, where a is the standard deviation of the yields in a given run range, are considered 
outliers. Two such runs were found out of 457 considered. The integrated luminosities of 
these two runs are 14.3 nb~^ and 258.3 nb~^. Further investigations of online operational 
conditions during these runs revealed no obvious hardware or trigger malfunctions. Since 
the probability is 1% that a data subsample of 258.3 nb~^ out of a total sample of 39.7 
pb~^ would have a yield different by > 4cr, both runs are included in the baseline cross- 
section measurement. The measurement is repeated without the outlier runs included and 
a systematic uncertainty assigned from the difference in the measurements. We find the 
uncertainty on the total cross section to be less than 1%. 

D. Trigger Efficiency 

For our measurement of the Level 1 di-muon trigger efficiency, we used J/ events that 
were taken with a high-pr single- muon trigger. At Level 1, this trigger requires a muon 
with pt greater than 4.0 GeV/c. In Level 3, a, J/t/j is reconstructed using the triggered high- 
er muon and a second muon which is not required to pass the Level 1 requirements. This 
second muon is then used to measure the Level 1 single-muon efficiency. The denominator of 
the efficiency measurement is the number of J/ip reconstructed using the Level 3 track and 
muon information. These J/t/j candidates must have a mass between 2.7 and 3.6 GeV/c^, 
a di-muon opening angle of A0o < 130°, and a separation in zq of less than 5 cm between 
the candidate's tracks. The probe-muon track must have at least 20 COT axial-layer hits 
and 16 COT stereo-layer hits, a CMU r-0 match of x^(Ar0) < 9, and a track |-2o| < 90 cm. 



Tracks are excluded if they pass within 1.5 cm of the center of any of the COT wire planes 
in any of the axial layers in order to avoid the inefficient region caused by wire supports. For 
the probe muon to pass the Level 1 trigger, the associated Level 3 track must be matched 
to an XFT track and the Level 3 CMU stub must be matched to a Level 1 CMU stub that 
lies within XTRP window. The resulting Level 1 muon-finding efficiency is shown in Fig. 7. 
The distribution is fit to the following function: 
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FIG. 7: The Level 1 CMU trigger efficiency as a function of muon pT- Points with error bars are 
measurement points. The solid line is the fitting result using the function described in the text. 
The dashed Unes are range used to determine the uncertainty. 



where freq is the normal frequency function: 



freq(a;) = 



1 



e~2* dt, 



(6) 



(7) 



E is the plateau efficiency, A is associated with the pr at which the efficiency is half the 
peak value, and R is the effective Gaussian resolution. We find E — 0.977 ± 0.002, A — 
1.1 ± 0.1 ( GeV/c)-\ and R = 0.28 ± 0.06 ( GeV/c)-\ 



To determine the uncertainty in the Level 1 trigger efficiency, while also taking into 
account the data fluctuations around the central fit as shown in Fig. 7, the range of the 
uppermost and lowermost fluctuations supported by the data are computed as follows: 
x'{pt) = X ± {\x — x\ + la) where x is the data value, x is the value returned by the 
fit and a is the uncertainty on the data. The x'ipr) distribution is refit using the function 
in Equation 6. The results are shown as dashed hues in Fig. 7. The di-muon Level 1 trigger 
efficiency is calculated on an event-by-event basis to take into account n-fx correlations. For 
each J/ip candidate, the Level 1 J/ip reconstruction efficiency is given by: 

eUApJ'') - eM') ■ eUpn (8) 

where eLi(p^) is the single muon Level 1 trigger efficiency given by Equation 6, and pf^'^ are 
the transverse momenta of the two muon candidates. The trigger's exclusion of pairs with 
nearby stubs is included as part of the geometric acceptance. The mean of the Level 1 di- 
muon trigger efficiency distribution in each J/ip transverse momentum bin is listed in Table 
L The maximum difference from varying the trigger efficiencies by one standard deviation 
independently for the two muons is listed as the uncertainty on the di-muon trigger in Table 
I. We find that the variation is within ±1.5% in all bins. 

The Level 3 reconstruction efficiency is dominated by the difference between the online 
and offiine tracking efficiency. A fast tracking algorithm is used for pattern recognition in 
the COT in Level 3. In the offiine reconstruction a more accurate tracking algorithm is 
combined with the result of the Level 3 algorithm to give a higher overall COT tracking 
efficiency. The Level 3 single-muon reconstruction efficiency as measured versus the offiine 
reconstruction algorithm is found to be constant for pr(A*) > 1.5GeV/c and is 

^L3/offline = ^'^^^ ± O.OOl(stat) ± 0.002{syst). (9) 

In the Level 3 trigger, the muons are required to be separated in zq by less than 5 cm. 
The efficiency eAzo of this cut is measured using J/ip candidates reconstructed in single- 
muon-trigger data samples where a Level 3 di-muon trigger was not required to acquire the 



data. The numbers of events that passed the Zo-separation criterion in the mass signal and 
sideband regions are examined. The cut is found to be 100% efficient with an uncertainty 
of 0.1%. The uncertainty is driven by the statistical limitations of the small data samples 
obtained from the single-muon triggers. 

E. Reconstruction Efficiencies 

The COT tracking efficiency was measured using a Monte Carlo track embedding tech- 
nique. Hits from simulated muon tracks are embedded into CDF Run II di-muon events. 
The distance resolution and hit- merging distance are adjusted so the embedded track has 
residuals and hit distributions matched to muon tracks in J/%1) data events. The efficiency 
of COT track reconstruction in di-muon events is found to be 

ecoT(PT > 1-5 GeV/c) = 0.9961 ± QmQ2{stat)to^^{{syst). (10) 

The absolute offline reconstruction efficiency of muons including stub reconstruction and 
matching stubs to tracks is measured using J/%Ij events from single-muon trigger samples 
where the J/%1) invariant mass is reconstructed from a triggered, fully-reconstructed muon 
and a second track. Tracks from the di-muon-mass signal region are projected to the muon 
chambers, and the efficiency of finding a matched stub is measured. For muons in the CMU 
fiducial region with Pt(a*) > 1.5GeV/c, the offline reconstruction efficiency is found to be 
independent of pr and is measured to be: 

e^j^u = 0.986 ± 0.003 ± 0.010. (11) 

To select clean CMU muons, the track-stub matching in the r-cf) plane is required to have 
X^(Ar0) < 9. The efficiency of this cut is found to have a weak dependence on p^: 

e^z = (1.0018 ± 0.0003) - (0.0024 ± 0.0001)p^. (12) 

The efficiency of the track-stub matching criterion (x^(Ar0) < 9) as a function of J/V' 
transverse momentum, obtained using an event-by-event weighting is listed in Table I. The 



TABLE II: Summary of J/ij; reconstruction efficiencies. 



J/ip Selection 


Efflciency 


Level 3 muon reconstruction 


6^3 = 0.997 ± 0.001 ± 0.002 


COT offline tracking 


ecoT = 0.9961 ± 0.0002;°;°°^^ 


Muon offline reconstruction 


ecMU = 0.986 ± 0.003 ± 0.010 


Muon zq position less than ±90 cm 


e^^ = 0.9943 ± 0.0016 


Di-muon zq separation less than 5 cm 


eAzo = 1-0 ± 0.001 


Total reconstruction e^ec 


~ ^L3 ' ^COT ' ^CMU ' • CAzo = 9^-^ ^ 2.7% 



systematic uncertainty on the weighted- average matching-cut efficiency is obtained by vary- 
ing the normahzation and slope in Equation 12 by one standard deviation. The change in 
the weighted average efficiency in each J/il) transverse momentum bin is found to be < 0.2%. 

Since the two muons originate from a common decay point, the efficiency of the track 
Zq cut is fully correlated for the two muons and is counted only once. The combined px 
independent COT-tracking, muon and Level 3 reconstruction efficiencies for J/%1) mesons is 
calculated to be 

Table II summarizes the pr-independent reconstruction efficiencies and those of the various 
muon selection cuts. 

V. J/ijj CROSS SECTION 

An event-by-event weighting is used to determine the J/ip yield in each Pt bin. Each 
event is weighted using the Level 1 single muon efficiency eLi(p^) and the efficiency of the 
track-stub matching criterion e^2(p^) applied to each of the two muons. The event is then 
corrected for the acceptance A{p^J^ -ly"^^^). The weight of each candidate event is given by: 

1/wi = tu{p^) ■ eLiiPr) ■ 

e^.{f^).eAv'^).A{p'4\y"'^)- (14) 



We fit the invariant mass distributions of tfie weigfited events, using tlie same sfiapes for 
signal and background as shown in Fig. 2, 3, 4 and 5. The number of signal events in each 
transverse momentum bin is determined from the area under the signal mass peak. The 
error on the corrected yield from the mass template fit, A'"(pt) corrected is given by: 



S{N{p' 



'T J corrected J 



\ i=0 



where Ng is the raw number of signal events in each momentum bin before weighting. In 
a similar fashion, the di-muon px distribution in each bin is weighted. The weighed pt 
distribution of the mass sideband subtracted events in the J/ip mass signal region is used 
to determine the mean pt value for each transverse momentum bin. 
The J/ijj differential cross section is then calculated as follows: 

■DtJll \ ^(Pr) corrected ' - A!) 

— .BHJIi,^,,) = ^—j^^^^^, (16) 

where da/dpr is the average cross section of inclusive J/ip in that pt bin integrated over 
I y{J/ip) \< 0.6, A' is the correction factor for y smearing defined by Equation 4, erec is 
the combined Level 3 and offline tracking and muon reconstruction cfflciency, J Cdt is the 
integrated luminosity, and Ap^ is the size of the px bin. 

The cross-section values obtained with statistical and pr-dependent uncertainties are 
listed in Table 111. 

An uncertainty of +0.1% on the momentum scale is extracted by comparing the re- 
constructed J/V' mass as shown in Fig. 1 to the world averaged value of 3.09688 ± 
0.00004 GeV/c^ [22]. The 3 MeV/c^ difference is attributed to an underestimation of the 
energy loss in the silicon detector due to an incomplete accounting of the material at the 
time the data sample used in this analysis was processed. The +0.1% uncertainty from 
the momentum scale corresponds to an uncertainty on the differential cross section as 
d{da / dpx) / dpT x 0.1%. Using the values in Table III, the first derivative of the differ- 
ential cross section is calculated and the momentum scale uncertainty on the cross section 
in each bin estimated. The effect was found to be small, the largest negative deviation being 
—0.08% and the largest positive deviation being +0.7%. 



TABLE III: The differential J/if) cross section times tlie brancfiing fraction Br = Br{J/'tp — /jl/i) 
as a function of pr for \y{J/iJ;)\ < 0.6. For each measurement, the first uncertainty is statistical and 
the second uncertainty is systematic. The systematic uncertainties shown are the pT dependent 
uncertainties only. The fully correlated pr independent systematic uncertainty in each bin is 7%. 



Pt{J/'^) ( GeV/c) Mean pr Mean p^, ^ ■ Br (nb/( GeV/c)) ^ • Br (nb/( GeV/c^)) 



0.0 — 

0.25 - 


0.25 

- 0.5 


0.15 

0.39 


0.027 

0.16 


9.13 ± 0.6(stat)lj!,'7(syst) 
28.1 it 1.53j^ g 


36.5 ± 2A{stat)_2 g{syst) 
37.4 ± 2.0^2 2 


0.5 — 


0.75 


0.64 


0.42 


45.3 ± 1.912 ° 


36.2 ± 1.5l^'g 


0.75 - 


- 1.0 


0.89 


0.79 


59.3 ± 2.0_2 9 


33.9 ± l.l^^ g 


1.0 — 
1.25 - 


1.25 
- 1.5 


1.13 
1.38 


1.29 
1.91 


69. D ± l-9_3 2 

73.4 ± 1.7_3 g 


31.0 ± 0.8_]^ g 

26.7 ± 0.61j 3 


1.5 — 


1.75 


1.63 


2.66 


75.2 ± 1.6t3'| 


23.2 ± 0.51J'q 


1.75 - 


- 2.0 


1.87 


3.52 


72.9 ± 1.4t| 3 


19.4 ± 0.4_g g 


2.0 — 
2.25 - 


2.25 
- 2.5 


2.13 
2.38 


4.53 
5.65 


69.1 ± 0.8^2 g 
67.3 ± 1.0_2 8 


16.3 ± 0.2+g * 
14.2 ± 0.2_g g 


2.5 — 


2.75 


2.62 


6.89 


57.6 ± 0.9 ± 2.6 


11.0 ± O.i ± 0.5 


2.75 - 


- 3.0 


2.87 


8.26 


52.0 ± 0.8 ± 2.4 


9.04 ± 0.13 ± 0.41 


3.0 — 


3.25 


3.12 


9.76 


43.6 ± 0.7 ± 1.9 


6.97 ± 0.10 ± 0.31 


3.25 - 


- 3.5 


3.38 


11.4 


37.3 ± 0.6 ± 1.6 


5.53 ± 0.08 ± 0.24 


3.5- 


3.75 


3.62 


13.1 


31.5 ± 0.5 ± 1.3 


4.34 ± 0.07 ± 0.18 


3.75 - 


-4.0 


3.87 


15.0 


26.2 ±0.4 ± 1.2 


3.38 ±0.05 ±0.15 


4.0- 


4.25 


4.12 


17.0 


22.5 ±0.4 ±1.0 


2.72 ±0.05 ±0.12 


4.25 - 


-4.5 


4.38 


19.2 


18.7 ±0.3 ±0.8 


2.13 ±0.04 ±0.09 


4.5- 


4.75 


4.62 


21.4 


16.1 ±0.3 ±0.7 


1.74 ± 0.03 ± 0.08 


4.75- 


-5.0 


4.88 


23.8 


13.3 ±0.3 ±0.6 


1.37 ± 0.03 ± 0.06 


5.0 - 


- 5.5 


5.24 


27.5 


10.3 ±0.15 ±0.42 


0.984 ±0.014 ±0.040 


5.5 - 


- 6.0 


5.74 


33.0 


7.28 ±0.12 ±0.29 


0.633 ± 0.010 ± 0.025 


6.0- 


- 6.5 


6.24 


38.9 


5.11 ±0.09 ±0.20 


0.408 ± 0.0069 ± 0.016 


6.5 - 


- 7.0 


6.74 


45.5 


3.54 ±0.07 ±0.14 


0.262 ± 0.0052 ± 0.010 


7.0- 


- 8.0 


7.45 


55.7 


2.27 ±0.03 ±0.10 


0.151 ± 0.0019 ± 0.006 


8.0- 


- 9.0 


8.46 


71.6 


1.14 ±0.02 ±0.05 


0.0668 ± 0.0011 ± 0.0028 


9.0- 


10.0 


9.46 


89.5 


0.622 ±0.013 ±0.025 


0.0327 ± 0.0007 ± 0.0013 


10.0- 


- 12.0 


10.8 


118 


0.278 ± 0.006 ± 0.011 


0.0126 ± 0.0003 ± 0.0005 


12.0- 


- 14.0 


12.8 


165 


0.103 ±0.003 ±0.004 


0.00398 ± 0.00013 ± 0.00015 


14.0 - 


- 17.0 


15.2 


233 


0.037 ± 0.002 ± 0.002 


0.00120 ± 0.00005 ± 0.00006 


17.0- 


- 20.0 


18.3 


336 


0.014 ±0.001 ±0.001 


0.00037 ± 0.00004 ± 0.00002 



Table IV summarizes the different contributions to the systematic errors applied to the 
cross-section measurement from acceptance calculations using a Monte Carlo simulation, 
the mass line shapes used to determine the yield, the trigger and reconstruction efficiencies. 



and the luminosity measurement. 



TABLE IV: Summary of systematic uncertainties in the inclusive J/ip cross-section measurement. 
The pT dependent uncertainties are listed in the top section of the table. In general, the pr depen- 
dent uncertainties increase with decreasing pT- The total is calculated from the pT independent 
sources only. 





Source 




Size 


Acceptance 


J/^ spin alignment 


±(2 


^ 5)% {pt) 


Acceptance 


Pt spectrum 


±(0 


^ 5)% {pt) 


Acceptance 


Detector material 


±(0.4 


^ 5)% {pt) 


Yield 


Mass fits 


(-1.3 


+9)% {pt) 


Yield 


Momentum scale 


(-0.1 ^ +0.7)% ipr) 


Luminosity 


CLC 




±6.0% 


Reconstruction 


Table II 




±2.8% 


Acceptance 


CMU simulation 




±1.0% 


Yield 


Data quality 




±1.0% 


LI trigger efficiency 


Table I 




±1.5% 


Total 


±6.9% 


®5{pT 


) 



The differential cross-section results with systematic and statistical uncertainties are dis- 
played in Fig. 8. The invariant cross section, da/dp^-Br{J/ip —>■ nn), with systematic errors 
is shown in Fig. 9. The results are also listed in Table III. 

We integrate the differential cross section to find the total J/^ production cross section: 

a{pp J/^X, I y{J/ij) \< 0.6) • Br{J/ip /x/x) 

= 240 ± l{stat)t^^{syst) nb. (17) 

The pr-dependent systematic uncertainties are summed and then added in quadrature with 
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FIG. 8: Inclusive J/V' cross section, da/dpr ■ Br{J/ip — ^ fi/i), as a function of J/V' Pt integrated 

over the rapidity range \y\ < 0.6. The differential cross section with systematic and statistical 
uncertainties added is plotted. This includes correlated uncertainties. 
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FIG. 9: The invariant J/ip cross section, da/dp^ ■ Br{J/ip — /x/x), as a function of p^{J/il)). The 
differential cross section is plotted. This includes correlated uncertainties. 



the fully-correlated uncertainty of 6. 



xstat 
tot (T rmBr 



70. 



A^bin 



= ^11 nb, (19) 

where A^bins is the total number of pr bins, is the statistical uncertainty in the cross- 
section measurement in the i*'* bin, 5^^^^(pt) is the systematic uncertainty on the measure- 
ment in each pt bin independent of the correlated systematic uncertainty of 6.9%, and © 
denotes addition in quadrature. After correcting for the Br{J/ip /i/i) = 5.88±0.10% [22], 
we find 

a(pp^J/^ljX,\y(J/jlj)\<0.6) 

= 4.08 ± 0.02{stat)toitisyst) /^b. (20) 

To compare with prior measurements where only the portion of the cross section for 
Pt{J/iP) exceeding 5GeV/c was measured [1], we also measure the integrated cross section 
of inclusive J/t/j with pt > 5GeV/c and \r]\ < 0.6 at ^/s — 1960 GeV. We find the cross 
section is 

a{pp J/^X,pTiJ/ip) > 5.0 GeV/c, | ri{J/ip) \< 0.6) • Br{J/^ iijj) 

= 16.3 ± Q.l{stat)t\^^{syst) nb. (21) 

We discuss the comparison of this result with earlier data in Section VII. 



VI. Hb J/ip FRACTION AND THE &-HADRON CROSS SECTION 

In general, the inclusive J/ip cross section contains contributions from various sources: 
prompt production of charmonium; decays of excited charmonium states such as ip{2S), Xci 
and Xc2', and decays of 6-hadrons. The charmonium states decay immediately. In contrast. 



6-hadrons have long lifetimes that are on the order of picoseconds. This implies that J/ip 
events from the decays of 6-hadrons are likely to be displaced from the beamline. We exploit 
this feature to separate J/ ^0 of decay products of 6-hadrons from that of prompt charmonium 
in the pt bins used in the inclusive J/ if) cross section calculation. 

To measure the fraction of events that are from displaced decay vertices, we use the 
subset of the J/%1) sample that includes those events for which both muon tracks from the 
J/V' satisfy high quality COT-SVX II track requirements. The track extrapolation from the 
path formed by the trajectory in the COT into the SVX II is described in Section II C 1. The 
total number of hits expected in the five layers of the SVX II is determined from the number 
of functioning and powered sihcon sensors intersected by the COT muon track. Tracks 
missing more than one expected hit in the SVX II are rejected. Both tracks are required to 
have a hit in the innermost layer of the SVX II and a hit in the second layer if the sensor 
intersected by the COT track is functioning. Corrections for energy loss in the SVX II are 
applied to the candidate muons based on a GEANT simulation of the material. The two 
muon tracks are constrained to come from a common space point. The probability of 
this 3-dimensional vertex fit is required to exceed 0.001. We find that 139200 ± 500 events, 
or about half of the total data sample, pass these criteria. While the data sample 
is reduced by the SVX II requirements, the momentum, angle, and vertex resolutions are 
substantially improved. 

The primary vertex, taken as the beam position in the r-(j) plane, is assumed as the 
point where 6-hadrons are produced. It is calculated on a run-by-run basis from a data 
sample taken using the inclusive jet trigger which has negligible contributions from charm 
and bottom decays so the beamline position can be calculated with no bias from detached 
decay vertices. The resolution of the primary vertex in the r-0 plane is limited by the 
~ 30 jJLva RMS spread in the size of the beam envelope. 



A. Measurement of the Fraction of J jip Events from 6-hadrons 



The J/^ from the decay of H\, J /ipX is hkely to be displaced from the primary vertex 
where 6-hadrons are assumed to be produced. The signed projection of the flight distance of 
J/'0 on its transverse momentum, Lxy, is a good measurement of the displaced vertex and 
can be used as a variable to separate J/ip of the iJ^ decay products from that of prompt 
decays. This method works well for events with high J/ip pt where the flight direction aligns 
well with that of the 6-hadron. For events with very low J/ip pr, the non- negligible amount 
of J/ip with large opening angle between its flight direction and that of the 6-hadron will 
impair the separation ability. Monte Carlo simulation shows that a reliable 6-fraction can 
be extracted using this method for events with J/t/j pt greater than 1.25 GeV/c. 

The Lxv is calculated as 



where L is the vector from the primary vertex to the J/ip decay vertex in the r-0 plane 
and priJ/ip) is the transverse momentum vector. To reduce the dependence on the J/ip 
transverse momentum bin size and placement, a new variable x, called pseudo-proper decay 
time, is used instead of Lxy, 



where the M{J /ip) is taken as the known J mass [22]. A Monte Carlo simulation is needed 
to model the distribution of x{J/i/j) from 6-hadron events. The Monte Carlo templates of the 
X distributions X^c(x,p^^) are generated for all J/ijj transverse momentum ranges and are 
directly convoluted with the value of the x resolution function measured in the data without 
allowing any of the parameters governing the shape of the Monte Carlo distributions to vary. 

1. The Likelihood Function 

An unbinned maximum likelihood fit is used to extract the 6-fraction, fs, from the data. 
The J/ip pseudo-proper decay time x, its error a, and the mass of the di-muon pair m^y^^ are 



L^y{j/i;)^L-pT{j/^)/\PT{jm 



(22) 



x^L^y{J/^)-M{J/^)/pT{J/i^), 



(23) 



the input variables. A simultaneous mass and lifetime fit is performed using a log-likelihood 
function (ln>C) given by: 

N 

ln£ = ^lnj^(x,m^^), (24) 

i=l 

where N is the total number of events in the mass range 2.85 < m^^ < 3.35GeV/c^. 

The mass and pseudo-proper decay time distribution is described by the following func- 
tion, 

J^{x,mi^i^) = /sig X J^sig(a;) x A4sig("^MM) 

+ (1 - /sig) X J^Bkg{x) X AlBkg(m^^), (25) 

where /sig is the fraction of signal J/ip events in the mass region, Tsig and Tskg are the 
functional forms describing the J/i/j pseudo-proper decay time distribution for the signal and 
background events respectively, and A^sig and A^Bkg are the functional forms describing the 
invariant mass distributions for the signal and background events respectively. We now 
describe these components of the likehhood fit in more detail. 

The function for modeling the J/ ip pseudo-proper decay time signal distribution consists 
of two parts, the Hj, — > J/ipX decay and prompt decay functions labeled Tb{x) and J^p{x) 
respectively: 

J'sigix) = [fB ■ :Fb{x) + (1 - Jb) ■ :Fp{x)] , (26) 

where is the fraction of J/ij: mesons originating in 6-hadron decays. We use the x 
distributions of accepted events from a Monte Carlo simulation as templates for the x 
distribution of 6-hadron events in data. The generated distributions are convoluted with a 
resolution function R{x' — x, sa) such that the Hf, — > J/ipX signal shape is given by 

J^b{x) = R{x - X, sa) ® X^^{x'), (27) 

where s is an overall error scale factor which represents the possible errors in determining the 
lifetime resolution and ® denotes a convolution. Prompt J/ if) mesons are produced at the 
primary vertex, therefore their observed displacement is described only by the resolution 



function Tp = R{x,sa). We find that R{x' — x,sa) is best described by a sum of two 
Gaussian distributions centered at x = 0. 

The background requires a more comphcated parameterization to obtain a good fit to the 
data outside the J/ip signal region. The pseudo- proper decay time background function is 
composed of four parts: the zero hfetime component, a positive slope exponential function, 
a negative slope exponential function, and a symmetric exponential function with both 
positive and negative slopes. The positive slope exponential function is chosen to model the 
background from other long lived 6-hadron events that produce opposite sign muons such 
as 6 — > cii~PX, c — > /I'^uX. The zero lifetime component is chosen to be the same shape 
as the resolution function. The symmetric and negative slope exponential functions are 
added to parameterize the remaining components of the background pseudo-proper decay 
time distributions which are from unknown sources. The background exponential tails are 
also convoluted with the resolution function. 

The background functional form is parameterized as follows: 



where /-t. sym is the fraction of the background distribution in the positive, negative and 
symmetric exponential tails respectively, A±_sym are the corresponding exponential slopes, 
and 9{x) is the step function defined as 9{x) = 1 for x > and 9{x) = for x < 0. It 
should be kept in mind that the background strongly depends on and m^^, and that the 
likelihood function incorporates a global fit over the full mass window shown in Fig. 2 to 
Fig. 5, including the J/ip peak and mass sidebands. 

The mass resolution used in the likelihood fit is better than that shown in Figs. 3-5 



^Bkg{x) = (!-/+-/_- fsyrajRix, Sa) 





(28) 



because of the addition of SVX II hits to the tracks. For the hkehhood fit, the di-muon 
mass shape Alsig is chosen to be simply the sum of two Gaussian distributions. The means 
of the Gaussian distributions are allowed to float independently: 

Msigirrif,^) = Gi{m^^- M^gm) 

+ /2 • ^2^^ - (M + D),r2aM). (29) 

The mass fit parameters are the mean M of the mass distribution, the width um of the 
first Gaussian distribution, the fraction /2 of the second Gaussian distribution, the shift D 
in the mean of the second Gaussian distribution, and the ratio r2 of the widths of the two 
Gaussian distributions. The mass background is modeled using a linear distribution. This 
fit is adequate for the SVX II constrained di-muon mass. The function used, normalized to 
unity over the mass range m"*'" to 777,"*"^ is: 

1 



"Vm mm 



+ M,iope(m,, - ^ ), (30) 

where Mgiope is the slope of the mass background distribution. The only fit parameter is 

-^slope- 

2. The Fits and Systematic Uncertainties 

The fits to the J/ip pseudo-proper decay time in three sample px ranges are shown in 
Fig. 10, Fig. 11 and Fig. 12. These data correspond to a subset of the data in the mass plots 
shown in Fig. 2 to Fig. 5 which satisfy the SVX II tracking requirements. The values of the 
6-fractions from the fits with statistical and systematic uncertainties for events with J/ip 
transverse momenta of 1.25 GeV/ c to 20.0 GeV/ c are listed in Table V, and the distribution is 
shown in Fig. 13. This measurement of the 6-fraction is used in Section VI B, in conjunction 
with the measurement of the inclusive J/ip cross section, to calculate the inclusive 6-hadron 
cross section. 
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FIG. 10: Fits to the J/ip pseudo-proper decay time in the range 1.25 < pr(/^Ai) < 1.5GeV/c to 
extract the fraction of events from long-lived 6-hadron decays. The solid line is the fit to all the 
events in the mass window of 2.85 to 3.35 GeV/c^, the dashed line is the fit to all signal events, the 
solid histogram is the fit to the portion of the signal events that are from 6-hadron decays and the 
dot-dashed line is the fit to background events including events in the invariant mass sidebands. 

The uncertainties on the 6-fractions are summarized in Table VI. In the table, percentage 
errors on the absolute value of 6-fraction are listed. Now we discuss the estimation of 
systematic uncertainties on the 6-fraction in detail. 

We have performed various tests to assess the accuracy of the likelihood procedure. The 
fit shapes for signal and background are histogramed into bins and compared to the binned 
data distributions. A Kolmogorov-Smirnov (K-S) test [39] is used to compare the fit and 
data distributions to estimate the quality of the fit. The distribution of K-S probability 
values for each fit in the different transverse momentum ranges is compared to the K-S 
probability distributions in a sample of Monte Carlo experiments. The distributions are 
found to be consistent. In addition to the K-S tests, the normalized residual, defined as 
the difference between the data and fit projection in the unit of one standard deviation of 
statistical error, is compared in every transverse momentum range. Firstly, the data and 
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FIG. 11: Fits to the J/ip pseudo-proper decay time in the range 5.0 < pxifJ-fJ') < 5.5GeV/c to 
extract the fraction of events from long-hved 6-hadron decays. The sohd hne is the fit to all the 
events in the mass window, the dashed line is the fit to ah signal events, the sohd histogram is the 
fit to the portion of the signal events that are from 6-hadron decays and the dot-dashed line is the 
fit to background events. 

fit projections are histogrammed using an unequal pseudo-proper decay time bin size so 
that the number of data events in each bin is more than 20 events to reduce statistical 
fluctuation. Secondly, the normalized residual distributions are examined. The means and 
widths of the distributions in all transverse momentum ranges are examined. We flnd no 
obvious discrepancies between the flt projection and data distributions. 

Monte Carlo samples are also used to determine the potential bias on the 6-fraction from 
the fltting procedure. The pseudo-proper decay time distributions and the invariant mass 
distributions from signal and background are used to generate a set of 500 statistically 
independent samples for each of the four pt bins of 1.25-1.5 GeV/c, 2.0-2.25 GeV/c, 5.0- 
5.5GeV/c and 10.0-12.0 GeV/c. Five different values of the 6-fraction, 5% to 13%, are 
assumed and the number of events in each pt bin is chosen to match the data. The fltted 
values of 6-fractions are found to agree with the generated values within 2% over the whole 
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FIG. 12: Fits to the J/^p pseudo-proper decay time in the range 12.0 < pxifJ-fJ-) < 14.0GeV/c to 
extract the fraction of events from long-hved 5-hadron decays. The sohd hne is the fit to all the 
events in the mass window, the dashed line is the fit to ah signal events, the sohd histogram is the 
fit to the portion of the signal events that are from 6-hadron decays and the dot-dashed line is the 
fit to background events. 

Pt bins. Thus the systematic uncertainties on the 6-fraction measurements due to fit bias 
are found to be less than 2%. 

The resolution function for the pseudo-proper decay time, R{x' — x, sa), is modeled by 
a double Gaussian function where the dominant Gaussian width is allowed to float and is 
determined by the fit to the data in each Pt{J/4') bin. Other parameters in the function 
are fixed to the values obtained from a binned fit to Lxy/cr(Lxy) averaged over all pxiJ/'ip)- 
The double Gaussian resolution function is not an exact description of the resolution func- 
tion shape but only an approximate parameterization of many different resolution effects. 
Therefore, to estimate the systematic uncertainty due to the resolution function modeling, 
the maximum range of values for the ratios of areas and widths of the two Gaussians sup- 
ported by the data are estimated. We find that the ratios of the second Gaussian to the 
dominant Gaussian vary from 1.5 to 2.5 in width and 0.05 to 0.15 in area. The systematic 
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FIG. 13: Fraction of J/V' from 6-hadron decays in the inclusive J/V' events of Run-II data as a 

function of J/V' transverse momentum. Error bars include both statistical and systematic errors. 

uncertainty on the 6-fraction from this source is largest in the lowest momentum bin, where 
the percentage error is as large as 8%, and decreases with increasing prl-Z/V')- 

In the Jjil) pseudo-proper decay time signal region, events are observed in the distribution 
at long positive and negative hfetimes that are not well described by the double Gaussian 
description of the resolution function. The source of these long lived "tails" is unknown. To 
estimate the systematic uncertainty on the long lived tails not modeled by the prompt signal 
double Gaussian, a box shaped function is added to the prompt J/ip x signal distribution 
in the range -2500 to 2500 yum. The height of the box is fixed in the fit using the number 
of events in the data that are in excess of the fit projection. The 6-fraction values returned 
from the fit with the box function are used to estimate the systematic uncertainty from 
the tails that are not modeled properly. We find the 6-fraction values decrease by about 
5% in the lowest momentum bins when the box shape is added to the prompt J/%1) distri- 
bution. The excess modeled by the box can also be assigned to the 6-hadron signal which 
causes a systematic increase. The change in the 6-fraction decreases in the higher transverse 



TABLE V: The fraction of J /if) events from decays of 6-hadrons and the corresponding acceptance. 
The first uncertainty on the 6-fraction is the statistical uncertainty from the unbinned hkehhood 
fit and the second uncertainty is the combined systematic uncertainties on the measurement of the 
6-fraction. The uncertainty on the acceptance is the combined statistical uncertainty from Monte 
Carlo statistics and the systematic uncertainty on the acceptance measurement. 
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momentum bins. 

The fit was repeated with the background shape changed such that only a positive 
and negative exponential is used with no symmetric exponential. The differences in the 
6-fractions observed are negligible. The background parameters are extracted from a fit to 
the sideband data distributions only, where the sidebands are chosen such that no significant 



TABLE VI: Sources of systematic uncertainties on the measurement of the 6-hadron fraction in in- 
clusive J/ if) decays as percentages of the absolute value. In general, the pT dependent uncertainties 
increase with decreasing px- 



Source 


Systematic 


uncertainty 


Resolution function model 


±(0.5 
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-2)% 
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-2)% 
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-7)% 
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±(0.5 


-3)% 
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-4)% 


Total 
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13)% 



contribution is expected from the radiative J/ if) tail. The fit is repeated in each bin with 
the values of the background parameters fixed to the values obtained from the sideband fit. 
No significant difference between the value of the parameters extracted is observed. The 
difference in the 6-fraction extracted using the parameters obtained from the sideband fit is 
taken as a systematic uncertainty. In the lowest and highest momentum bin the percentage 
difference on 6- fraction value extracted is 2-3%. The differences are less than 1% in all other 
bins. 

To study the dependence of the 6-fractions on the modeling of the 6-hadron spectrum 
used in the Monte Carlo, a fiat distribution in pt and y of the h production spectrum 
is used to regenerate the x distributions and the fits were repeated. The differences in 
the value of the 6-fractions extracted from the direct fit to Monte Carlo templates of x 
produced from an input spectrum that is uniform in prih) and y{h) are examined. The 
variation in the 6-fractions extracted in the range 1.25 to 2.0GeV/c are the largest, the 
maximum variation being an increase of 18% in the bin 1.5 to 1.75GeV/c. The uniform 
input spectrum is unrealistic, therefore the systematic uncertainty is taken as one-half of 
the size of the variation observed in the 6-fraction when the fiat model is used. We assign 



systematic uncertainties of 7%, 3% and 2% for measurements in the transverse momentum 
ranges of 1.25-3.0 GeV/c, 3.0-8.0 GeV/c and 8.0-20.0 GeV/c respectively. 

In addition, we examine the change in the 6-fraction extracted when varying the — > 
J/iIjX decay momentum spectrum while keeping the Hb production momentum spectrum 
fixed. Two decay spectrums, Hi, J/'^(direct)X and Hi, J/-^ (inclusive) X [40], are used 
for this. The percentage difference is found to be 2-4% in the lowest momentum bins and 
< 1% ioipT^J/i^) > 2.5GeV/c. 

The mix of hadrons and their respective lifetimes is a contributing factor to the shape of 
the J/ip pseudo-proper decay time distributions. To assess the systematic uncertainty due to 
the uncertainty on the 6-hadron average lifetime, we vary the average lifetime in the Monte 
Carlo by 11 /xm which is the size of the systematic uncertainty on the average 6-hadron 
hfetime measured at CDF in Run II. We find that the measured 6-fraction decreases in all 
transverse momentum bins when the lifetime is increased. The fractional decrease is 4% in 
the lowest momentum bin and less than 1% for bins with Pt{J/ 'ip) > 12 GeV/ c. The variation 
in the 6-fraction observed when the average 6-hadron lifetime is varied by ±11 /im is taken as 
a systematic uncertainty on the 6-fraction measurement. Table VI summarizes the sources 
of systematic uncertainties on the measurement of the 6-hadron fraction as percentages of 
the absolute values. 

B. Measurement of the Inclusive 6-hadron Cross Section 

Since J/ ip mesons from decays of bottom hadrons have a different average spin alignment 
than an inclusive sample of J/%1) mesons, we need to apply an acceptance correction to 
account for this difference. In previous CDF measurements, the effective value of the spin 
alignment parameter cteff of J/il) from 6-hadron decays was measured to be otesipT^J > 
4.0 GeV/c) = —0.09 ± 0.10 [5], where a^e is obtained by fitting cos^j/y,, the angle between 
the muon direction in the J/V' rest frame and the J/V' direction in the lab frame, to the 
functional form 1 -|- cteff • cos^ More recent measurement on the spin alignment was done 



using B — > J/il)X events collected at the T(4S') resonance. The BaBar experiment measured 
ae = -0.196 ± 0.044 for p* < 1.1 GeV/c and ae = -0.592 ± 0.032 for p* > 1.1 GeV/c [41]. 
Here the decay angle of the J/t/j is measured in the T(45') rest frame and p* is the total J/ip 
momentum measured in the T(4S') rest frame. 

We opt to use the more precise result from the BaBar experiment in the acceptance 
calculations for H^, — > J/t/jX events assuming it is apphcable to the CDF environment 
where 6-hadrons are produced in fragmentation with a large momentum range instead of 
produced at a fixed momentum as in T decays [15]. First, Monte Carlo events are generated 
to have the J/ip helicity angle distributions in the 6-hadron rest frame predicted from as 
values according to their p* values. Then, values of the spin alignment parameter agff for 
events in each J/ ip pt bin are obtained from fitting the cos 9j distributions of these Monte 
Carlo events. The systematic errors on are obtained by varying the input values of ckb 
in the process according to measurement errors. This process gives a result consistent with 
previous CDF measurement, albeit with smaller uncertainties. For example, a new and more 
precise value of ctefr = — 0.13±0.01 for the J/t/j events withpT(-'^/'0) > 4.0 GeV/c is obtained 
from this process. Finally, the acceptance values, as hsted in Table V, are calculated from 
the Monte Carlo events generated with the derived spin alignment parameters in each J/ip 
Pt bin. 

The differential 6-hadron cross sections are calculated in a similar way as that for the 
inclusive J/ip. The J/ip yields in each px bin listed in Table I are multiplied with the b- 
fractions to obtain the corresponding Hi, ^ J/V' yields. The new acceptance values hsted 
in Table V are used while the J/ i/j reconstruction efficiencies and luminosity value stay the 
same. Most of the systematic uncertainties in the inclusive J/ip cross-section calculation 
carry over here without change except for those from the J/ip spin alignment on the ac- 
ceptance which are estimated using errors on CKeff- In addition, the uncertainties from the 
6-fractions are also included in the systematic errors. The J/ip from 6-hadron inclusive 
cross-section results with statistical and systematic uncertainties are shown in Table VII. 
The differential cross section with all statistical and systematic errors added is plotted in 



Fig. 14. A recent QCD theoretical calculation using a fixed order (FO) calculation with re- 
summation of next-to-leading logs (NLL) [42] is overlaid. We discuss further the comparison 
with theoretical calculations in Section VII. 
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FIG. 14: Differential cross-section distribution of J/^ events from the decays of 6-hadrons as a 
function of J/^jj transverse momentum integrated over the rapidity range \y\ < 0.6. The crosses 
with error bars are the data with systematic and statistical uncertainties added including correlated 
uncertainties. The solid line is the central theoretical values using the FONLL calculations outlined 
in [42], the dashed line is the theoretical uncertainty. 

An integration of the differential 6-hadron cross-section results in Table VII gives the 
total cross section 



The systematic uncertainty quoted includes the fully correlated uncertainty of 6.9% obtained 
from the inclusive J/ip cross-section measurement. We correct the integrated cross section 



a{pp ^ H,X,pT{J/i^) > 1.25 GeV/c, \y{J/^)\ < 0.6) 



■Br{Ht, J/ipX) ■ Br{J/i/j ^j) 



lQA±Q.?,{stat)t\l{syst) nb. 



(31) 



TABLE VII: The inclusive J/ipX and prompt J/ip differential cross sections as a function 

of transverse momentum of the J/ip with statistical and pT dependent systematic uncertainties. 
The cross section in each pT bin is integrated over the rapidity range \y{J/ip)\ < 0.6 The fully 
correlated systematic uncertainty, systjc = 6.9%, from the measurement of the inclusive J/ip cross 
section needs to combined with the px dependent systematic uncertainties. 
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33.4 ±0.5 ±2.8 
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3.75 


-4.0 


3.87 


3.28 ±0.14 ±0.16 


22.9 ± 0.3 ± 1.6 


4.0- 


-4.25 


4.12 


2.45 ±0.11 ±0.15 


20.1 ±0.4 ± 1.5 


4.25 


-4.5 


4.38 


2.22 ±0.10 ±0.11 


16.5 ±0.3 ±1.2 


4.5- 


-4.75 


4.62 


1.99 ±0.09 ±0.10 


14.1 ±0.3 ±1.0 


4.75 


-5.0 


4.88 


1.84 ± 0.08 ± 0.10 


11.5 ±0.3 ±0.8 


5.0 - 


- 5.5 


5.24 


1.38 ±0.05 ±0.06 


8.92 ±0.13 ±0.52 


5.5 - 


- 6.0 


5.74 


1.07 ±0.04 ±0.05 


6.21 ±0.10 ±0.37 


6.0 • 


-6.5 


6.24 


0.817 ±0.031 ±0.038 


4.29 ±0.07 ±0.24 


6.5- 


-7.0 


6.74 


0.610 ± 0.025 ± 0.026 


2.93 ±0.06 ±0.17 


7.0- 


- 8.0 


7.45 


0.447 ±0.014 ±0.022 


1.82 ±0.02 ±0.11 


8.0 - 


- 9.0 


8.46 


0.246 ±0.009 ±0.010 


0.894 ±0.015 ±0.047 


9.0 - 


- 10.0 


9.46 


0.149 ±0.007 ±0.006 


0.473 ±0.010 ±0.024 


10.0 • 


- 12.0 


10.8 


0.074 ± 0.003 ± 0.003 


0.204 ± 0.004 ± 0.010 


12.0 ■ 


- 14.0 


12.8 


0.034 ± 0.002 ± 0.001 


0.069 ± 0.002 ± 0.003 


14.0 - 


- 17.0 


15.4 


0.0143 ± 0.0009 ± 0.0007 


0.023 ±0.001 ±0.001 


17.0 ■ 


- 20.0 


18.3 


0.0062 ± 0.0006 ± 0.0004 


0.0078 ± 0.0006 ± 0.0006 



extracted above for the branching fraction Br{J j'^ /x/i) = 5.88 ± 0.10% [22] to obtain 

(7(pp ^ E^, ^ J/i;,pT{J/i^) > 1.25 GeV/c, \yiJ/^)\ < 0.6) 

= 0.330 ± Q.QQ5{stat)t^Q^oll{syst) /xb. (32) 

We also extract the prompt J/ip cross section by subtracting the cross section of Hh 



J/ipX from the inclusive J/ip cross section. This calculation is applied to all J/ip with 
Pt > 1.25GeV/c where we are able to extract the 6- fraction. The results are shown in 
Table VII and in Fig. 15. The systematic uncertainties on the prompt J ftp cross section are 
taken to be the uncertainties on the inclusive cross section added in quadrature with the 
uncertainties on the measured 6-fractions. We find the integrated cross section of prompt 




Pt(J/v) GeV/c 



FIG. 15: The inclusive J/^l^ cross section as a function of J/^l^ pT integrated over the rapidity range 
\y\ < 0.6 is plotted as points with error bars where all uncertainties have been added. The hatched 
histogram indicates the contribution to the cross section from prompt charmonium production. 
The cross-hatched histogram is the contribution from decays of 6-hadrons. 

J/ip to be: 

aipp ^ J/iljpX,pTiJ/^p) > 1.25 GeV/c, |y(J/^)| < 0.6) 

= 2.86 ± OM{stat)toitisyst) /ib, (33) 

where J/ ipp denotes a prompt J/ ip and where we have corrected for the J/ip ^ /i/i branching 
fraction. 

The differential fe-hadron cross section as a function of pT{Hb) is extracted from the 
measured differential cross sections of Hb —>■ J/ipX by utilizing the decay kinematics of 



charmonium produced in 6-hadron decays. 

The procedure starts with the calculation of contributions to the cross section of 6-hadrons 
in a given Pt{Hi,) bin from J/t/j events in the range 1.25 < pt{J/iP) < 20GeV/c, where we 
measured the fractions of J/ip mesons from b decays. Since 6-hadrons with as little as zero 
momenta produce J/ip mesons with momenta as large as 2 GeV/c, the measured cross section 
in this pt{J/iP) range is sensitive to the complete pt{Hi,) spectrum. The total contribution 
to the cross section in the i*^ bin in pT{Hb) from events in the accessible Pt{J/'4^) range, 
labeled as the raw cross section ai{raw), is given by 

N 

ai{raw) = ^Wijaj{J/tjj), (34) 
j=i 

where aj{J/ip) is the cross section of J/il) mesons from H^, in the Pt{.J/4^) bin and Wij is 
the fraction of events in the i*'* pt^H^) bin relative to the total in the j*^ PriJ/'ip) bin. 
The sum of the weights Wij in each pt{J/iP) bin is normalized to 1. The raw cross section 
is corrected for the acceptance due to the limited J/ip Pt range to obtain the differential 
6-hadron cross section, ai{Hb), in the i^^ pT{Hb) bin, 

'Ji{Hb) = — = Ti > (35) 

where is the fraction of bottom hadrons in the i*'^ priHf,) bin that give rise to a J/V' with a 
transverse momentum in the range 1.25 to 20 GeV/c and rapidity in the range \y{J/^p) \ < 0.6. 
Monte Carlo simulations are used to calculate the weighting factors, Wij, and acceptance 
correction factors, /*. In the simulation, the decay spectrum of — > J/ipX obtained from 
references [40] and [41] is used. The calculation is repeated in an iteration process: at each 
pass the input production spectrum used in the Monte Carlo is the spectrum measured in 
the previous iteration and a comparison is made between the input and output spectrums. 
The process terminates when the comparison reach the precision limit. This procedure 
is found to be insensitive to the initial production spectrum shape. 
The statistical uncertainty in each pt{Hi,) bin is given by: 



SstaticTiiHb)) = 77 



N 



T.'^AtMJm- (36) 



The systematic uncertainties are taken as just the simple weighted sum of the systematic 
errors from the differential Hh —> J/ip cross-sections measurements, 

1 ^ 

Ssyst{(^i{Hb)) = -F-Yl ^iAyst{(TjiJ/ij))- (37) 
J'^ i=l 

The extracted differential cross section of 6-hadrons over the transverse momentum range 
from to 25GeV/c is shown in Fig. 16. The cross section has been corrected for the 
branching fractions, Br{Hb — >• J/ipX) = 1.16 ± 0.10% and Br{J/ip — > jjbjji) = 5.88 ± 
0.10% [22], and divided by two to obtain the single 6-hadron differential cross section. We 
integrate the differential cross section extracted above to obtain the single 6-hadron inclusive 
cross section. We find the total inclusive single 6-hadron cross section is 

a{pp HhX, \y\ < 0.6) = 17.6 ± QA{stat)tli{syst) iih. (38) 
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FIG. 16: Differential cross-section distribution of 6-hadron production as a function of 6-hadron 
transverse momenta. The crosses with error bars are the data with systematic and statistical 
uncertainties added, including correlated uncertainties. The solid line is the central theoretical 
values using the FONLL calculations outlined in [42] , the dashed line is the theoretical uncertainty. 



VII. DISCUSSION 



We have measured the inclusive J/ip and 6-hadron cross sections in pp interactions at 
y/s — 1960 GeV in the central rapidity region of \y\ < 0.6. For the first time, the cross 
section has been measured over the full transverse momentum range (0-20 GeV/c). 

For comparison to Run I measurements at y/s = 1800 GeV [1], we consider the cross- 
section measurements in the range pt('//'0) > 5.0 GeV/ c and pseudo-rapidity \r]{J/ip)\ < 0.6. 
We measure the inclusive J/t/j cross section at ^/s — 1960 GeV to be 

a{pp J/il)X)ig6o ■ Br{J/ip jjLjj) 

= 16.3 ± Q.l{stat)t\^^{syst) nb. (39) 

The CDF Run I measurement at = 1800 GeV was found to be 

a{pp J/ipX)isoo ■ Br{J/ip nn) 

= 17.4 ± Q.l{stat)tli{syst) nb. (40) 

We measure the cross section of J/ip events from Hi, decays with Pt{J/iP) > 5 GeV/c and 
\v{J/i^)\ < 0.6 to be 

a{pp ^ HbX)igeo ■ Br{Hb ^ J/i^X) ■ Br{J/ij; ^ nn) 

= 2.75 ± OM{stat) ± 0.2Q{syst) nb. (41) 

The equivalent Run I measurement [1] was found to be 

a{pp ^ HbX)i8oo ■ Br(Hb J/^X) ■ BriJ/i) i^ii) 

= 3.23 ± 0.05{stat)toil{syst) nb. (42) 

Although the Run II J/ il: and 6-hadron cross sections are measured at a higher center-of-mass 
energy, and it is expected that the production cross sections increase by approximately 10%, 
the Run I and Run II measurements are consistent within measurement uncertainties. The 
ratio of the Run II to Run I differential 6-hadron cross-section measurements as a function 
of Pt{J/iP) is shown in Fig. 17. No difference in the shape of the cross section is observed. 
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FIG. 17: Ratio of the differential cross-section distributions of J/V' events from the decays of 
6-hadrons as a function of J/ij) transverse momentum from CDF Run I and Run II. 

In Fig. 18, the B'^ differential cross section previously measured by CDF at — 
1800 GeV for \y\ < 1.0 [9] is compared with our newer measurement of the inclusive b- 
hadron differential cross section at ^/s — 1960 GeV extracted from the measurement of the 
cross section of J/i/j events from 6-hadron decays. For the purpose of this comparison, the 
CDF Run II inclusive 6-hadron cross section is multiplied by the fragmentation fraction of 
mesons, where the result from LEP experiments is used [43]. In addition, the Run II 
6-hadron inclusive cross section is scaled up by a factor of 1.67 to extend the measurement to 
\y\ < 1.0 where we have assumed the rapidity distribution is uniform in the region \y\ < 1.0. 
As shown in Fig. 18, we find good agreement between the Run II extracted measurement of 
the 6-hadron cross section and the direct measurement of the B'^ cross section in Run I. 

In Fig. 14 and Fig. 16, we compare our measurement to a QCD calculation of the 6-hadron 
cross section by Cacciari et al. [42]. This calculation uses a fixed-order approach with a 
next-to-leading-log resummation and a new technique to extract the 6-hadron fragmenta- 
tion function from LEP data [20, 42]. The single 6-hadron cross section from this FONLL 
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FIG. 18: The differential cross-section distributions of 6-hadrons as a function of 6-hadron mo- 
mentum from the measurement of the meson cross section in CDF Run I [9] and 6-hadron 
inclusive cross section extracted in this analysis (Run II). The differential cross sections showed 
are integrated over the rapidity range \y{Hh)\ < 1.0. 

calculation using the CTEQ6M parton distribution functions [44] is = 16.8^5;o /ub 

which is in good agreement with our measurement of 17.6 ± 0A(stat)^2'l{'Wst) /^b. 

We also compare this result to the QCD calculation described in reference [19]. This cal- 
culation employs a factorization scheme where the mass of the quark is considered negligible 
and a different treatment of the 6-hadron fragmentation function is used. The cross-section 
calculation in [19] is repeated using y/s = 1960GeV/c and the MRST2001 parton distri- 
bution functions [46]. The central value of the calculated cross section integrated over the 
rapidity range \y\ < 0.6 and priJ/i^) > 5.0GeV/c is a^pp H^X, \y\ < 0.6) • Br{Hh — *■ 
J/ipX) ■ Br{J/il) fifi) = 3.2 nb [45] which is in good agreement with our result of 
3.06 ± 0.04(stat) ± 0.22{syst) nb 

A more complete discussion of the changes in QCD calculations can be found in refer- 
ences [18, 20, 42]. Updated determinations of proton parton densities and bottom quark 



fragmentation functions have brought the QCD calculations into better agreement with the 
CDF measurements of the total 6-hadron cross section and the 6-hadron px distribution. 

VIII. SUMMARY 

We have measured the inclusive central J/V' cross section in pp interactions at ^/s ~ 
1960 GeV. The cross section has been measured over the full transverse momentum range, 
for the first time. We find the integrated inclusive J jil) cross section in the central rapidity 
range to be 

a{pp^Jli)X,\y{Jli,) |<0.6) 
= 4.08 ± 0.02(stoi)12J^(sysi) //b, (43) 

after correcting for Br{J jii nn) = 5.88 ± 0.10% [22]. 

Using the long lifetime of 6-hadrons to separate that portion of the J/ ■0 cross section that 
is from decays of 6-hadrons, we have measured the cross section of J/ ip mesons from 6-hadron 
decays for J/ip transverse momenta greater than 1.25GeV/c. The integrated Hj, J/ipX 
cross section, including both hadron and anti-hadron states, is 

a(pp ^ H,, H, ^ J/'ijjX,pT(J/ip) > 1.25 GeV/c, |y(J/V')| < 0.6) 

= 0.330 ± Q.QQ5{stat)to°oll{syst) ^h, (44) 

after correcting for the branching fraction Br{J/ip — > /x/x) = 5.88 ± 0.10% [22]. 

The measurement of the J/ ip cross section from 6-hadron decays probes 6-hadron trans- 
verse momenta down to zero. We have extracted the first measurement of the total central 
6-hadron cross section in pp collisions from the measurement of the 6-hadron cross section 
with J/t/j transverse momenta greater than 1.25GeV/c using Monte Carlo models. We find 
the total single &-hadron cross section integrated over all transverse momenta to be 

a{pp HbX, \y\ < 0.6) = 17.6 ± OA{stat)tli{syst) //b. (45) 
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